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ABSTRACT 

Low resolution SDSS/SEGUE spectra have been used to study the behavior of 
the strengths of the CN and CH molecular bands in stars at different evolutionary 
stages of the open cluster NGC 6791. We find a significant spread in the strengths 
of the CN bands, more than twice that expected from the uncertainties, although 
the bimodalities observed in globular clusters are not clearly observed here. This 
behavior, is observed not only among red clump objects but also in unevolved 
stars such as those in the main sequence and lower red giant branch. In contrast, 
not all the stars studied show significant scatter in their CH strengths. 

Subject headings: stars: abundances — open clusters and associations: individ- 
ual: NGC 6791 



1. Introduction 



From the pioneering work of lOsbornl (ll97lf ). variations in the strength of CN molecular 
bands have been obs erved in globular cluster stars of similar luminosities a t different evo- 



lutionary stages (e.g. ICannon et al.lll998t ICohenlll999l : lHarbeck et al.ll2003l ). The strength 



of these molecular bands depends on temperature, surface gravity, and chemical compo- 
sition. Therefore, the fact that stars of similar luminosity and at the same evolutionary 
stage have different CN band strengths denotes a different chemical composition. These 
variations i n CN band strengt h are often anti-corr elated with the CH band intensity at 
4300A (e.g. ICannon et al.lll998t lHarbeck et al.ll2003 ). Moreover, detailed analyses based on 



high-resolution spectroscopy reveal that variations in CN strengths are accompanied by the 
existence of anticorrelations of the chemical abundances of other light element s (Na, O, Al, 
and Mg), but not for heavier elements such as Fe or Ca dCarretta et al]l2009al lb[). 
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Several hypotheses concerning intrinsic a nd extrinsic factors have been proposed to ex- 
plain these results (see iPancino et al .11201 Obi , for a recent review). However, the detection 



of mu ltiple evolutionary sequences in almost all clusters properly studied (e.g. iPiotto et al. 
20071 ) indicat es that the most promising one is the so-called self-enrichment scenario. First 
proposed by iDa Costa fc Cottrelll (119801 ) . this hypothesis assumes the existence of two or 
more subsequent stellar generations. The first stars formed pollute the gas from which 
the second population is formed. It seems that almo st all globular cluste rs are massive 
enough to host at least two stellar populations (e.g. ICarretta et al.l 120101 ) . It is there- 
fore important to investigate what happens in less massive systems such as open clus- 
ters, which are usually younger, less massive and more metal-rich than globulars. How- 
ever, stu dies based on open clusters using both the strength of CN bands in low-resolution 



spectra ( INorris fc Smith! Il985l ; iHufnagel et al.l |l995; Ma rtell fc Smith! 120091) and chemical 



abundances derived from high- resolution spect ra (Ide Silva et al. 



2009 



Smiljanic et al.l 12009 



Pancino et al.ll2010ar ICarrera fc Pancindl201ll ) have not found the same trends observed in 
globulars. 

There is an intriguing stellar system, NGC 6791, with a mass intermediate between 
those of open and globular clusters. NGC 6791 is an ideal target for exploring the existence 



2006 



Carraro et al.l 120061) . and high metal content 



Platais et al. 


2011 


), old age ( 


~+0.4 dex 


Carraro et al. 


2006 



8 Gyr; e.g. 



Carretta et al.ll2007f ) differentiate this system from other open clusters. Furthermore, 
it is located about 1 kpc above the Galactic plane and the ecce ntricity of its orbit is greater 
than the typical value for old open clusters (IBedin et al.ll2006|). Previous determ inations of 



NGC 6791's orbit are compatible with an extragalactic origin (jCarraro et al.ll2006l ). However, 
proper motions derived from Hubble Space Telescope observations suggest that this cluster 
was formed near the Galactic bulge with a higher initial mass ( IBedin et al.l 120061 ). The 
total mass of NGC 6791 has decreased over time owing to repeated crossings of the denser 
parts of the disk. The color-magnitude diagram of NGC 6791 shows an unexpectedly 
wide red giant branch (RGB) and main-sequence (MS) in the region of the turn-off. Some 
authors ha ve explained these features in terms of the existence of differential line-of-sight 
reddening ( IPlatais et al.l 120111). An alternative explanation, however, is the existence of an 
age spread of about 1 Gyr ( jTwarog et al.l 120111 ). As well as its red clump (RC), typical of 
a metal-rich population, NGC 679 1 also has a handful of extreme blue ho ri zontal branch 



stars whose origin is still unclear ( IPlatais et al.l l201ll ; iBuzzoni et al.l 120121 ) . iKalirai et al. 



(120071) associate these wi th strong winds due to the high-metallicity at the end of the RGB, 



while IBedin et al.l ( 120081 ) suggest that they are the consequence of a high binary fraction. 
Recent investigations of globular clusters associate these extreme blue horizontal branch 
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stars with the exis tence o f a He -rich population (e.g. iGratton et al.ll2010l ) although this issue 
remains unsolved. iJaned ( 19841 ). using DDO photometr y, found no sign s of va riations in the 
compositions of the RGB stars analyzed. In contrast, iHufnagel et al.l (119951 ) found strong 
indications of intrinsic differences among the CN molecular band strengths in the RC stars 
of NGC 6791. Investigations of the composition of the stars in NGC 6791 b ased on high 



resolution spectroscopy ha ve found no evidence of chemical inhomogeneities ( lOriglia et al. 



20061 ; ICarretta et al.l 120071). altho u gh an unexpected sc atter in the Na abundance has been 
reported by ICarretta et al. ( 2007 ). Geisler et al. ( 2012 ) were recently the first to report the 
existence of anticorrelations in the Na and O abundances in NGC 6791, similar to those 
observed in globular cluster stars. The goal of this paper is to investigate the behavior of 
the CN and CH molecular band strengths in the stars of NGC 6791 at different evolutionary 
stages. 

This paper is organized as follows: the observational material is described in Section [2] 
molecular indexes are defined in Section [3j the CN and CH distributions are obtained and 
analyzed in Section HJ and the main results of this work are discussed and summarized in 
Section O and Section El respectively. 



2. Observational Material 



In the frame work of the Sloan E xtension for Galactic Understanding and Exploration 
(SEGUE) survey (lYanny et al.l 120091 ) . low resolution (R ~ 2000) spectra in the wavelength 
range of 3600-9200A were obtained for more than 1000 stars observed in two plates in the 
line of sight of NGC 6791 . They were observed to validate th e SEGUE Stellar Param- 
eter Pipeline (SSPP, lAllende Prieto et al.l 120081 ; iLee et al.ll2008al Jbl). The reduced spectra, 
together with the atmosphere parameters, radial velocities, etc., derived for observed stars 



were made public in the eighth Sloan Digital Sky Survey (SDSS) data release (lAihara et al. 



20111 ). The procedure for obtaining this information can be summarized as follows. Raw 
sp ectra are first reduced by the SDSS spectroscopic reduction pipeline, described in detail 
by lStoughton et al.l (120021 ). which provides flux- and wavelength-calibrated spectra, together 
with initial determinations of the radial velocities and spectral types. More accurate ra- 
dial velocities are calculated in a subsequent step by SSPP, together with determinations of 
metallicity, effective temperature, and surface gravity. Two examples of the spectra used in 
this paper are shown in Figure [TJ 



To select members of NGC 6791 from all the observed stars we followed a similar proce- 
dure as that described by lSmolinski et al.l ( I2011al ). We first reject those stars located farther 
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than 5^50 from the cluster centei@ (Fig. [2]). This radius was chosen to ensure that stars 
selected were among those observed in the two photometric samples used throughout this 
paper and described below. We then discard those stars whose radial velociti es (Fig. are 



not w ithin the range (K)ngc679i ± lc 5 assuming (V^.)ngc679i=-46±10 km s 1 ( ICarrera et al. 



20071 ). In the same way, we reject those stars whose metallicities (Fig. [3]) do not satisfy 
[Fe/H]±2cr with [Fe/H] = +0.3±0.15 dex (ICarrera fc Pancinoll201ll . and references therein). 
A total of 104 stars were selecte d as NGC 6791 member s according to these criteria. Al- 
though we used the same data as ISmolinski et al.l feOllaf ). our sample is bigger because we 
use a larger radius to select the cluster members. 

The SDSS data release also provides informat ion about the magnitudes of the tar- 



gets obtained with th e SDSS photometric pipeline (ILupton et al.ll2002l ; iTucker et al.l 12006 



Davenport et al.l 120071 ). However, this pipeline was not designed to handle the high-density 
crowded fields typical of the central areas of cl usters. Photome try from SDSS images of stellar 
clusters, i ncluding NGC 6791 was obtained bv kn et al.l koO$ ) using the DAOPHOT/ALL FRAME 
packages ( ]Stetsonlll987l . Il994j ) . The same packages have been used by lStetson et al.l (120031 ) to 
derive high-quality broadband BVI photometry for this cluster from 1764 CCD images. The 
resulting color-magnitude diagrams for NGC 6791 in each photometry system are shown in 
the left and right panels of Figure HI respecti vely. Owing to the hi gher quality and lower un- 
certainty of the BVI photometry derived by lStetson et al.l (120031 ). we use these magnitudes 
in this paper. The position of stars selected in these diagrams have been used to separate 
them into four groups according to their evolutionary stages: 71 MS (filled circles); 14 lower 
RGB (1RGB, crosses); 6 upper RGB (uRGB, filled triangles); and 6 RC (open stars). Five 
stars were discarded because they do not fall within any of these four groups (open cir- 
cles). The position of the RGB bump in an 8 Gyr old and Z ~ 0.04 isochrone (gray dashed 
line), selected from the BaSTI library (jPietrinferni et al.ll2004j ). was used to separate stars in 
the 1RGB and uRGB. We divided RGB stars in these two subgroups, because no pollution 
from material synthesized in the stellar interior is expected before the stars reach the bump. 
Coordinates, identification numbers, and magnitudes for selected stars are summarized in 
Table HJ 



1 This value is 1.25 x r/i and 0.25 x r t where =4.'4 and r t =23.', respectively (jPlatais et al.ll201li) . 
2 a 2 ooo = 19 h 20 m 53 s <5 20 oo = +37°46f3 
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Index definition 



The star-to-star light element abundance variations within globular clusters is investi- 
gated from the variations in the strength of the CN absorption bands around 3839 and 4142 
A, respectively, and the CH band around 4300 A (e.g. lNorris fc Freemanlll979uHarbeck et al. 
20031 : iKavser et al"1l2008[ IPancino et aDhoiObh . The strength of each molecular band is mea- 
sured with a spectral index defined as the magnitude of the difference between the integrated 
flux within a wavelength window containing the given feature and the integrated flux inside 
a nearby window, or windows, used to define the continuum. Several definitions of these 
indexes, used to measure the strength of each molecular band as a function of the luminosity 
class of the targets, can been found in the literature. This is particularly true for the 3839 
CN band, where the presence of other temperature-dependent absorption lines, such as the 
at 3839 A, complicates the definition of the continuum window. Since our targets cover 
a wide luminosity range, we used two differe nt indexes to measu re the strength of the CN 
band at 3839 A. The first one was defined by iNorris et al.l ( )198ll ) to sample RGB stars as: 



53839 



N 



-2.5 log 



3846-3883 



3883-3916 



The second one was defined by lHarbeck et al.l (120031 ) to sample properly MS objects as: 



S3839 H = -2.5 log' 



3861-3884 



3894-3910 



(2) 



The 4142 CN band becomes useful for samplin g CN variations in the cas e o f the relatively 



metal-r ich ([Fe/H]~-0.7 dex) globular cluster 47 ( INorris &: Freeman! 1 19 791 ). IPancino et al. 



( 12 01 Obi ) show that this band is less sensitive to CN variations for metal-poor MS stars. 
However, since NGC 6791 is much more metal-rich than globular cluste rs, we also studied 
the b ehavior of the the strength of this molecular band determined as (INorris Sz Freeman 
1979h : 



54142 = -2.5 log- 



4120-4216 



4216-4290 



(3) 



Finally, the G CH band at ~4300 A is typically used for sampling the carbon abundances 
in comparison with the CN bands which are used as i ndicators o f N abundances. We measure 
the strength of this band using the index defined by iLed ( 119991 ) : 
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Ci/4300 = -2.5 log 



F, 



4270-4320 



0.5^4230-4260 + 0.5^4390-4420 



(4) 



The adopted windows for each index ha ve been overp l otted i n Figure [TJ For each 
index, the uncertainty has been calculated as in lPancino et al.l (J2010bj) assuming pure photon 
(Poisson) noise statistics in the flux measurements. The indexes and errors determined 
together with photometric magnitudes for each selected star are listed in Table [TJ The 
median uncertainties for each index and evolutionary stage are listed in Table |2j 

As mentioned above, two instrumental configurations were employed for the NGC 6791 
region. A total of 33 stars were observed in both configurations although neither of them met 
the criteria used above to select NGC 6791 members. In any case, they are useful for checking 
the homogeneity of our data and provided an additional estimate of the uncertainties. To do 
this, each index was measured separately in the two spectra for each star. The medians of the 
differences between the values obtained for each index for all stars are: A(S , 3839at) = 0.00 ± 
0.02; A(S3839 H ) = 0.01±0.02; A(S4142) = 0.000±0.006; and A(C#4300) = -0.006±0.005. 
These differences are lower than the uncertainties obtained from the Poisson statistics. 



Hufnagel et al.l ( 119951 ) studied the CN and CH molecular band strengths in several RGB 



and RC stars of NGC 6791 from low resoluti on spectra (R< 1 000). Our sample has eight 
stars in common with that study. Although iHufnagel et al.l (119951 ) did not use the same 
index definitions used here, a comparison between both sets of measures reveals systematic 
effects due to d ifferent instrum e nts us ed in each case. The values obtained here against those 
determined by IHufnagel et al.l ( 119951 ) for each index have been plotted in Figure [5j A clear 
linear correlation is observed for the CN indexes, particularly in the case of the S4142 index. 
This is because we u sed the same windows to define the continuum and band intensity as 
Hufnagel et al.l ( 119951 ). although they obtained their indexes in a slightly different way. For 
each star we have calculated the difference between the value obtained here and that obtained 



by IHufnagel et al.l (119951 ). We have computed the standard deviation of the differences of all 
stars for given indexes. The values obtained are: 0.039, 0.036, 0.010, and 0.015 for S3839h, 
53839^, 54142 and CH4300, respectively. These values are of the order of the median of 
the uncertainties of each index for RGB or RC stars (see Table [2]). 



4. CN and CH distributions 

The run of the four molecular indexes defined above as a function of V magnitude are 
shown in Figure [6J For MS stars, the strength decreases as the V magnitude increases. Since 
temperature increases with magnitude for MS objects, the observed trend is explained by 
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the more efficient formation of CN and CH molecules at lower temperatures. For 1RGB 
objects, the behaviour is more random: as we ascend the RGB towards brighter magnitudes 
for V >15 the indexes >S3839# and 54142 increases, the index 53839^ slightly decreases, 
and C7f4300 remains almost constant. From V ~15 a possible turnover is observed in the 
strengths of the indexes studied. We will come back to this point in Section |5j 

4.1. Main Sequence and lower Red Giant Branch 

First, we are going to focus on those stars in the MS and 1RGB. As was explained 
above, according to stellar evolution models, it is expected that the chemical composition of 
their atmospheres may reflect the initial conditions of the molecular gas cloud from which 
the stars were formed. Before using CN and CH bands as chemical composition indicators, 
the temperature and gravity dependence of molecular band strengths should be removed. 
Different approaches can be found in the literature to removing this dependence, such as 
fitting the lower envelope or the median ridge line of a given molecular index as a function 
of color or magnitude. However, due to the low number of objects sampled in the case of the 
1RGB, this step should be performed with caution. To surmount this problem we adopted the 
following procedure. A linear least-squares fit was performed on half of the stars (randomly 
selected) in the 1RGB. This procedure was repeated 10 3 times using different random subsets 
each time. The final zero-point and slope for each index were obtained as the median of the 
zero points and slopes calculated in each individual test. The same procedure has been used 
for MS stars. With this procedure we tried to minimize the uncertainties due to the low 
number of object studied and the influence of the points on the edges. The final linear fit 
adopted in each case is shown as dashed lines in Figure E] and are listed in Table |3J We 
calculated the corrected pseudo-indexes, denoted by a 5 preceding the corresponding index, 
as the difference between the index and the adopted linear fit in each case. 

The histograms obtained for each corrected index are shown in Figure [7| for the MS 
(bottom) and 1RGB (upper), respectively. In each case, generalized histograms (solid lines) 
have been derived by assuming that each star is represented by a Gaussian probability 
function centered on its corrected index value whose a is equal to the uncertainty in the 
determination of the index. There is no hint of bimodalities for any of the four indexes in 
the case of the MS. However, the distribution widths of the two S3839 indexes seem larger 
than that expected from the median uncertainties, as shown in the bottom of each panel, 
which in both cases is about 0.04 (Table [2]). In fact, if we fitted each generalized histogram 
with a single Gaussian (dashed lines) the a values obtained are 0.083±0.001 and 0.091±0.001 
for S3839h and 53839^, respectively. These values are twice as great as the uncertainties. 
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In contrast, narrow distributions are obtained for 5*4142 and CHA300. Although the median 
uncertainties are similar to those of S3839 indexes (~0.04). In this case, the a obtained from 
the fit of a single Gaussian are 0.046±0.001 and 0.045±0.001, respectively. The residuals 
between the generalized histogram (solid lines) and the single Gaussian fitted (dashed lines) 
have been plotted in inset panels. They suggest that the distributions are not well reproduced 
by a single Gaussian. 

For 1RGB objects, the S4142 and CH4300 indexes, with a median uncertainty of ~0.03, 
also show a narrow distribution. They are reasonably reproduced by single Gaussians (dashed 
lines) with a ~0.037±0.001 and 0.031±0.001 for 54142 and CF4300, respectively. However, 
the residuals suggest that the generalized histograms have a wider dispersion than that 
obtained from the single Gaussian fit. In the case of two S3839 indexes a bimodality is sensed 
both in the normal and generalized histograms. Again, we have tried to fit each generalized 
distribution with a single Gaussian. In this case, the a values obtained are 0.086±0.001 and 
0.080±0.001, for 53839// and 53839jv, respectively. These values are almost three times the 
median uncertainties (~0.03) in each case. Again, the residuals of the difference between 
the generalized histogram and the single Gaussian fit suggest that a single Gaussian does 
not properly reproduce the observed distributions. 

It can be argued that the wider distributions observed in the case of the 53839 indexes 
are due to our assumption of a linear dependence of the molecular band strength with tem- 
perature and gravity. This approximation may be wrong in the case of the MS turn-off and 
the base of the RGB. For this reason we have repeated our analysis eliminating 1RGB stars 
fainter than V~17.25 mag and MS objects brighter than V~17.75 mag. Moreover, we have 
rejected MS stars fainter than V~18.75 mag owing to the larger uncertainties in the index 
determinations. The distributions of each index for the stars at each evolutionary stage are 
shown in Figure El They have been obtained following the same procedure as described pre- 
viously but using the restricted samples to perform the least-squares fits in order to remove 
the temperature and gravity dependence. In this case, the signs of bimodalities observed 
in the histograms of 53839 indexes for 1RGB stars seem clearer. However, these signs are 
smoothed when the generalized histograms are obtained taking into account the uncertain- 
ties. As before, the distributions are wider than those expected from the uncertainties. 
However, in this case the same behavior is also observed in the case of 54142 and Cif4300 
indexes. We have fitted each generalized histogram with a single Gaussian (dashed lines) 
and plotted the residuals in the inset panels. Again, the residuals suggest that the distribu- 
tions are not well reproduced by a single Gaussian with the exception of the CH4300 index. 
Therefore, we conclude that the behaviors described above for each index are not influenced 
by the assumption of a linear dependence of the molecular index strengths with temperature 
and gravity for stars near the MS turn-off and in the base of the RGB. Moreover, it seems 
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that the strengths of the 54442 band have the same behavior. On contrary, the CH band 
shows no spread. 



4.2. Red Clump 



As discussed above, iHufnagel et al.l ( 119951 ) studied 31 RC stars in a similar way to that 
described here. They found inhomogeneties of the CN bands but no signs of bimodalities as 
observed among globular clusters. I n our sample, we have six stars in the RC region, five 
of them having also been studied by IHufnagel et al.l (119951 ). In the case of RC stars, the V 
magnitude is not useful for removing the temperature and gravity dependence since they all 
have similar values. For this reason, in left panels of Figure [9] each index has been plotted 
against the B — / color. A linear fit has been performed for each index in order to remove 
the temperature dependence. The very small number of stars observed in the RC do not 
allow us to follow the same procedure used in the case of MS and 1RGB stars. In each case, 
the corrected index has been obtained as the difference between the index and the value of 
the linear fit. The histograms of the distributions obtained for each corrected index have 
been plotted in the right panels of Figure [9j The two S3839 indexes seem to have a bimodal 
behavior. However, this bimodality is smoothed when the generalized histogram is obtained 
taking into account the uncertainties (the solid lines in the right panels of Fig. [H]). In any 
case, the generalized histograms seem to be wider than that expected from the uncertainties. 
To investigate this point, we fitted a single Gaussian to each generalized histogram (dashed 
lines). The a values of the single Gaussian fitted are 0.053±0.005 and 0.056±0.005 for 
S3839h and 53839^, respectively. Therefore, the distributions ob tained are twice as wi de 
as is expected from uncertainties. A similar result was obtained by IHufnagel et al.l (119951 ) . 



In the case of the S4142 and CH4300 indexes, the distributions obtained are well repro- 
duced by a single Gaussian. Although the S4142 distributions seems slightly wider (~0.03) 
than that expected from the uncertainties. On the contrary, the CH4300 distribution has a 
width (~0.02) similar to the uncertainty. Th e lack of variati o ns in the CH band strength for 
RC stars in NGC 6791 was also reported by IHufnagel et al.l (119951 ). 



5. Discussion 



Significant variations and bimodalities in the strengths of the CN molecular bands, which 
are always anticorrelated with the CH stren gths, have been widely reported in globular clus- 
ter stars at different evolutionary stages (e.g. lCannon et al.lll998l ; ICohenlll999l ; lHarbeck et al. 



- 10 - 



20031 ; iKayser et al.ll2008t iPancino et al.ll2010bt ISmolinski et al.l l2011bl ). This result is ex- 
plained by the existence of different star-to-star C and N abundances. The rec ent discovery 
of se veral evolutionary sequences in almost all globular clusters studied (e.g. iPiotto et al. 



20071 ) has associated these chemical composition inhomogeneties with the existence of sev- 
eral stellar populations in each system. 

In contrast, studies on open clusters, which are more metal-rich, yo unger, and less 
mass i ve than the ^lobular s , have not been obser ved similar trends (e.g. iNorris fc Smith 
19851 ; iHufnagel et al.l Il995t iMartell fc Smith! 120091 ) . This is interpreted in terms of open 
clusters being formed by a single stellar population. However, similar abundance variations 
to those seen among globular stars would produce smaller scatter in the CN band strengths 
of open cluster stars owing t o their higher metallicities, with solar or above solar metal 
content (IHufnagel et al.lll995l ). Therefore, the same trends observed in globulars would be 
more difficult to detect in open clusters. 

In this sense, NGC 679 1 is a key system beca use it has a mass intermediate between 
globular and open clusters. IHufnagel et al.l ( 119951 ) reported an unexpected scatter among 
the strengths of the CN bands in the RC stars of NGC 6791. We confirm this spread not 
only in evolved RC stars but also among MS and 1RGB objects. Therefore, the different 
abundances responsible for this spread should be present in the molecular gas cloud from 
which NGC 6791 stars were formed. For example, the scatter observed in the strengths of 
CN band at 3839 A is twice that expected from the uncertainties. Moreover, some signs of 
bimodalities are sensed in the distributions of the indexes used to measure the strength of 
this band among 1RGB and RC stars. However, this result should be treated with caution 
owing to the the small number of stars studied in each case. We do not detect the CH-CN 
band strength anticorrelations observed in globular clusters, and all the stars studied seem 
to have a very similar CH band strength. 

It is difficult to explain the CN strength dispersion observed in terms of m i xing e ffects 
since this trend is also observed among unevolved stars. Moreover, iGeisler et al.l (120121 ) have 
recently reported an intrinsic dispersion of Na and O abundances among NGC 6791 stars 
in the whole RGB. These stars follow the same Na-0 anticorrelation observed in globular 
clusters. For analogy with globular clusters, these results can be explained by the existence 
of multiple stellar populations, NGC 6791 being the first open cluster to show this feature. 
In fact,the unexpectedly wide RGB and MS observed in the NG C 6791 color-magnitude 



diagram indicates an age spread of about 1 Gyr (jTwarog et a 



could also be due to the existence of differential reddening (IPlatais et al.l 120 111 ). 



20111). However, t hese features 



Recent investigations have demonstrated that light element abundances are correlated 
with the colors of RGB stars if near ultraviolet filters (i.e. Johnson U or SDSS u) are used 
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in globular clusters (e.g. Marino et al.l 120081 : Milone et al.l boiol . 120121 : lLardo et all 120121 ) 



although the details of this behavior are still unclear. An accurate wide-field photometric 
study including U filter would therefore provide information about the radial distribution 
of the different stellar populations, which are ke y to investigating the process of the cluster 
formation an d early chemical enrichment (e.g. iDecressin et al.l 120071 ; iD'Ercole et al.l 120081 ; 



Renzinill2008f ). 



Finally, the decrease in the strengths of the CN bands with magnitude for V < 15 
needs further consideration although it is beyond the scope of thi s paper. A simi l ar tren d 
has been observed in some more metal-poor globular clusters by ISmolinski et al.l ( 1201 lbl ). 
This turnover may not be unusual since it has been observed by other tra cers of CN ban d 
strengths, s uch as DDO photom etry, in other ope n cluster; e.g . NGC 188 f lMcClurelll974j ). 
NGC 2682 J Janes fc Smith! ll984T ). and NGC 7789 Jjaneslll977[ ). However, this question has 
not been studied in detail by any of these authors. A detailed study of DDO colors in stars 
of relatively metal-rich globular clusters shows that this decrease in CN strength is observed 
only for CN-rich stars and not for CN-weak objects. They compared their observations 
with colors obtained from synthetic spectra and concluded that, although it is expected that 
the mixing with material synthesized in the interior by the CN cycle should increase the 
N abundances, and therefore the CN strengths, the depletion of C would be so large that 
it would produce the decrease in the CN strengths observed, even though C is the lesser 
contributor to the formation of the CN molecule. In any case, a detailed analysis of this 
issue is necessary in order to understand the process that produces it. 



6. Summary 

We have studied the strengths of the CN and CH bands at 3839, 4142, and 4300 A in 
stars at different evolutionary stages in the open cluster NGC 6791. This system is one 
of the most massive open clusters known, and according to its orbit, it should have been 
more massive at the time of its formation. For this analysis, we used low-resolution spectra 
(R ~ 2000) obtained in the framework of the SEGUE project within the SDSS. Our main 
results are as follows: 

• A significant spread in the strengths of the CN molecular band at 3839 A is observed 
not only in evolved red clump stars but also among objects in the main sequence and 
lower red giant branch. This dispersion is at least twice as great as that expected 
from the uncertainties. This result is obtained by the two indexes used to measure the 
strength of this band. This scatter is not observed in the strengths of the CN band at 
4142 A. 
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• The distributions of the strengths of the two CN bands studied are not well reproduced 
by a single Gaussian even for main sequence stars. In fact, signs of bimodalities appear 
when the obtained distribution is compared with a single Gaussian. 

• No significant dispersion is observed in the strengths of the CH band at 4300 A. In 
fact, the distributions of the strengths of these indexes are relatively well reproduced 
by a single Gaussian. This is particularly true in the case of red clump stars. 

• The strengths of the three molecular bands studied decreases as the magnitude of the 
star is brighter in the upper part of the red giant branch above the red clump. 
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Fig. 1. — Example of two CN-weak (gray line) and CN-strong (black line ) spectra for two 
stars in the lower RGB. The IDs of these stars in the IStetson et al.l (120031 ) reference system 
are 14109 and 11885, respectively. The windows used for each index have been overplotted. 
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Fig. 2. — Stars w i th av ailable photometry obtained by I An et al.1 (120081 light gray points) 
and lStetson et all (120031 dark gray points) in the NGC 6791 region. Stars observed spectro- 
scopically by SEGUE and selected as NGC 6791 members have been plotted as black circles. 
Gray circles denote SEGUE stars discarded as cluster members. The dotted line denotes the 
tidal radius while the dashed line shows the radius used here to select NGC 6791 members 
(see text for details). 
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Fig. 3. — Radial velocity histograms of all stars observed by SEGUE (gray) and those located 
within 5f5 (black). Dashed lines denoted the radial velocity limits used. [Fe/H] distributions 
of stars selected from their radial velocity and from their [Fe/H] have been plotted in the 
inset panel with gray and black lines, respectively. 
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Fig. 4. — Color-magnitude diagrams from lAn et al.l ( 120081 . left) and IStetson et al.l ( 12003 
right). The dashed line is an 8 Gyr old isochrone from the BaSTI library (see text for 
details). Stars at different evolutionary stages have been plotted with different symbols: MS 
(filled circles); 1RGB (crosses), RC (open stars), and uRGB (filled triangles). Open circles 
are stars selected as NGC 6791 members but not included in any of the evolutionary stages 
defined. 
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Fig. 5. — Comparis on between the indexes measured here with those calculated by 
Hufnagel et al.l (119951 ) for stars in common. 
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Fig. 6. — Run of the four molecular indexes versus V magnitudes. Stars at different evo- 
lutionary stages have been plotted with different symbols as in Figure HI Dashed lines are 
the linear least-squares fit to those stars in the MS and 1RGB separately used to correct 
the temperature and gravity dependence (see text for details). Uncertainties obtained as 
the median of the errors of each index determination for each evolutionary stage have been 
plotted in the bottom part of each panel. 
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Fig. 7. — Histograms and generalized histograms (solid lines) for each index obtained for 
1RGB (top) and MS (bottom) stars, respectively. As a comparison, each generalized distri- 
bution has been fitted with a single Gaussian (dashed lines) and the residual between them 
are shown in inset panels. The residuals suggest that a single Gaussian does not properly 
reproduce the observed distributions. 




Fig. 8. — As Figure [7] but for 1RGB stars with V<17.25 and MS objects in the range 
18.75<V<17.75. 



-24 - 



-0.28 
-0.3 
-0.32 

-0.34 

0.06 
0.04 
0.02 


-0.02 
0.7 

0.65 
0.6 

0.55 
0.5 

0.35 
0.3 

0.25 
0.2 

0.15 



I | I I I I | I I I I | I I I I | I I I I | I I I I 



ir J" 



tr 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 



tr 



1 ir' " 



1 1 rT i 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 



ir 



tr X 



tr x 



- tr 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 



T 

tr -- 



I i i i i I i i i i I i i i i I i i i i I 



CH4300 



i i i i I i i i i I i i i i I i i i i I i i i i 
0.01 


-0.01 

-o.aoo.2 



A 



i ■ i ■ i 



■. i . i . i .- 



i i i i I i i i i I i i i i I i i i i I i i i i 



S4142 




S3839 



H 




0.2 


-0.2 




-0.2 0.2 



I i i i i I i i i i I i i i i I 



2.6 2.65 2.7 2.75 2.8 2.85 



-0.2 -0.1 0.1 0.2 0.3 



B-I 



Fig. 9. — Left: run of the strength of the molecular indexes of RC stars versus their B — I 
colors. A line has been fitted (dashed line) in order to remove the temperature and gravity 
dependence in order to obtain the corrected pseudo-indexes. Right: Normal (histogram) 
and generalized (solid line) distributions of each corrected pseudo-index. Dashed lines are 
the best single Gaussians fitted to each generalized distribution. The residuals between the 
single Gaussian and the generalized distribution are shown in inset panels. 



Table 1. Index measurements for the sample stars. 



Plate 


Fiber 


RA 


Dec 


g 


g-r 




V 


V-I 


S3839_ff 


S3839jv 


S4142 


Cff4300 


Region 


2800 


150 


19:21:04.28 


+37:47:18.84 


14.70 


1.42 


11814 


13.85 


1.66 


0.01+0.02 


0.33+0.02 


-0.07+0.02 


-0.30+0.02 


uRGB 


2800 


151 


19:21:14.52 


+37:46:32.81 


17.89 


0.94 


14109 


17.32 


1.14 


0.05+0.04 


0.51+0.05 


-0.15+0.04 


-0.32+0.04 


1RGB 


2800 


152 


19:21:06.70 


+37:48:08.31 


18.63 


0.71 


12451 


18.18 


0.94 


0.04+0.07 


0.29+0.07 


-0.26+0.06 


-0.32+0.06 


MS 


2800 


154 


19:21:01.45 


+37:48:05.11 


16.47 


1.14 


11006 


15.78 


1.31 


0.23+0.03 


0.55+0.03 


-0.05+0.02 


-0.31+0.02 


1RGB 



Note. — Table [T] is published in its entirety in the electronic edition of the Astrophysical Journal. A portion is shown here for guidance regarding its form 
and content. 
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Table 2. Median uncertainties for each index and evolutionary stage. 



Region 


S3839 H 




54142 


CH4300 


MS 


0.041±0.002 


0.045±0.002 


0.038±0.002 


0.038±0.002 


1RGB 


0.031±0.002 


0.034±0.003 


0.028±0.002 


0.027±0.002 


uRGB 


0.023±0.005 


0.025±0.001 


0.022±0.003 


0.021±0.002 


RC 


0.024±0.001 


0.026±0.001 


0.022±0.001 


0.020±0.001 



Table 3: Coefficients of the linear fit used to construct the corrected pseudo-indexes and 
marked as dashed lines in Figs M and EH respectively. 



Index 


Main-Sequence a 


lower Red Giant Branch a 


Red Clump b 




zp slope 


zp 


slope 


zp 


slope 


S3839 H 


0.23±0.06 -0.008±0.003 


0.67±0.20 


-0.03±0.01 


-1.39±0.25 


0.59±0.03 


S3839tv 


-0.46±0.07 0.049±0.004 


0.42±0.17 


0.01±0.01 


-0.94±0.26 


0.57±0.04 


S4142 


-0.64±0.02 -0.020±0.001 


1.08±0.06 


-0.070±0.004 


-0.78±0.03 


0.296±0.005 


CH4300 


-0.53±0.01 -0.011i0.001 


-0.32±0.04 


0.001±0.001 


-0.597±0.005 


0.108±0.001 



"Obtained as index = zp + slope x V. 
^Obtained as index = zp + slope x (B — I). 



